arXiv:1503.04972vl [math.CA] 17 Mar 2015 


Some new inequalities for the gamma function 

Xiaodong Cao 


Abstract 

In this paper, we present some new inequalities for the gamma function. The main tools 
are the multiple-correction method developed in [3 [7] and a generalized Mortici’s lemma. 


1 Introduction 


Duo to its importance in mathematics, the problem of finding new and sharp inequalities for 
the gamma function and, in particular for large values of x 



( 1 . 1 ) 


has attracted the attention of many researchers (see [SlElElEKIllOliailKITKIH] and the 


references therein). Let’s recall some of the classical results. Maybe one of the most well-known 
formula for approximation the gamma function is the Stirling’s formula 



See, e.g. [1, p. 253]. The following two formulas give slightly better estimates than Stirling’s 
formula, 



(Burnside [5], 1917 ), 


(1.3) 


(Gosper [TO], 1978). 


(1.4) 
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Ramanujian [22] proposed the claim (without proof) for the gamma function 


(1.5) 


r(x + 1 ) = 0F + 4x2 + X + 


where 0^; —>■ 1 as x —>■ +oo and ^ < 0^; < 1. This open problem was solved by Karatsuba[T3]. 
Thus (11.51) provides a more accurate estimate for the gamma function (see Sec. 2 below). 

In this paper, we will continue the previous works laizi, and introduce a class of new 
approximations to improve these inequalities. 

Throughout the paper, the notation \k(A:;x) denotes a polynomial of degree A; in x with all 
coefficients non-negative, which may be different at each occurrence. Let {an)n>i and {bn)n>o 
be two sequences of real numbers with 7 ^ 0 for all n € N . The generalized continued fraction 


T = bo + 


ai 


bi + 


= bo + 


ai 02 

^ 1 + b2+ 


• • • = 60 K 


—0 


b2 + 


is defined as the limit of the nth approximant 


— = bo + is. — 

Bn k—1 bf^ 


as n tends to infinity. See [ 2 , p.l05]. 


2 A generalized Mortici’s lemma 


Mortici |14] established a very useful tool for measuring the rate of convergence, which says that 
a sequence {xn)n>i converging to zero is the fastest possible when the difference (x^ — Xn+i)n>i 
is the fastest possible. Since then, Mortici’s lemma has been effectively applied in many paper 
such as igiTllIZllIB]. The following lemma is a generalization of Mortici’s lemma. 

Lemma 1. If liuix^+oo f (x) = 0, and there exists the limit 


(2.1) lim x^ (/(x) —/(x + 1)) = / G M, 

21—>-+00 

with A > 1, then there exists the limit 

( 2 . 2 ) lim x^~^f{x) = 

x—>-|-oo A — i 

Proof. It is not very difficult to prove that for x > 2 


(2.3) 


1 

(A — l)x'^“^ 




1 

(A-l)(x-l)^-i' 


2 









For e > 0, we assume that I — e < {f{x) — f{x + !))</ + £ for every real number x greater 

than or equal to the rank Xq > 0. By adding the inequalities of the form 

(/ - e)^ < f{x) - fix + 1) < (Z + 




m—1 


(2.4) 
we get 

m—1 m—1 

(2.5) il-e) Y, -^—^<fix)-fix + m) < {1 + e) Y + j)^ 

for every x > Xq and m > 1. By taking the limit as m ^ oo, then multiplying by x^~^, we 
obtain 


( 2 . 6 ) 


il-e)x^ ^Y 


1 


j=0 


(x+jy 


< x^-^fix) < il + e)x^-^Y 


1 


j=0 


(x+jy 


It follows from (j2.3p that 
(2.7) 


I — £ I + £ 


Now by taking the limit as x —)• +oo, this completes the proof of the lemma at once. □ 

An example Let’s consider the Ramanujan’s asymptotic formula (1.5). Let the error term 
Eix) be defined by the following relation 


( 2 . 8 ) 


r(x + 1) =-/tt + X + {l + E{x)). 


It follows readily from the recurrence formula r(x + 1) = xr(x) that 

(2.9) In (1 + E{x)) — In (1 + E{x + 1)) = — 1 + x In ^1 + 

1 8 (x + 1)^ + 4(x + 1)^ + (x + 1) + ^ 

6 8x^ + 4x2 + X + ^ 

By using the Mathematica software, we expand the right-hand function in the above formula as 
a power series in terms of 1 /x: 

(2.10) In (1 + Eix)) - In (1 + Eix + 1)) = + O(^). 

Thus, by Lemma 1 we have 

( 2 . 11 ) 

Noting that limu_,.o = 1, one get hnally 

( 2 . 12 ) lim x^Eix) = 

x->-i-cxD 11520 

Remark 1. Just as Motici’s lemma. Lemma 1 also provides a method for finding the limit of a 
function as x tends to infinity. 


liin x^ln(l+.F(x)) =-3^. 
x-)-i-oo 11520 
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3 Gosper-type inequalities 

In this section, we use an example to illustrate the idea of this paper. To this end, we introduce 
some class of correction function (MCfc(x))fc>o such that the relative error function Ek{x) has 
the fastest possible rate of convergence, which are defined by the relations 

(3.1) T{x + l) = ^/^{^ i + MCfc(a:) • exp(£’fc(x)). 

If lima;_>.+oo x^ f{x) = / / 0 with constant /i > 0, we say that the function f{x) is order and 
write the exponent of convergence fi = /i(/(x)). Clearly if fj.{Ek{x)) = fj-k, we have the following 
asymptotic formula 


(3.2) V{x + l) = V^ 



+ i + MCfc(x)- (l + 0(x-'^'=)). 


X +00. 


Let us briefly review a so-called multiple-correction method presented in our previous paper laE]. 
Actually, the multiple-correction method is a recursive algorithm, and one of its advantages is 
that by repeating correction process we always can accelerate the convergence, i.e. the sequence 
{fi{Ek{x)))k>o is a strictly increasing. The key step is to find a suitable structure of MCfc(x). In 
general, the correction function MCfc(x) is a hnite generalized continued fraction (see [7] or ()3.8h 
below) or a hyper-power series (see [6] or (14.7p below) in x. 

It is not difficult to see that (3.1) is equivalent to 

(3.3) lnr(x + 1) = ^ ln(27r) -b x (Inx - 1) -b ^ In (x -b MCfc(x)) -b Ek{x). 


By the recurrence formula r(x + 1) = xr(x), we have for x > 0 


(3.4) 


Ek{x) - Ek{x -b 1) 


— 1 -b X In 


1 

1 H— 

X 


1 j (x -b 1) + g + MCfc(x + 1) 

2 x + g-bMCfc(x) 


Now by taking the initial-correction function MCo(x) = and using Mathematica software, 
we expand Ek{x) — Ek{x -\-1 ) into a power series in terms of 1/x: 


(3.5) 


Eo{x) - Eo{x -b 1) 


— ^ kq 17 — 945/to — 810/to Aq 
x^ 540x^ 

-641 -b 33120/to - 12960/t2 -b 43200/toAo -b 25920/toAg /I 
12960x5 


The fastest possible function Eq{x) — Eq[x -b 1) is obtained when the hrst two coefficients in the 
above formula vanish. In this case, we find /to = ^, Ao = |g, and 


(3.6) 


Eo{x) - Eo{x -b 1) 


5929 

1166400x5 


+ 0 
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By Lemma 1, we can check that 


(3.7) 


lim x'^Eq{x) 

X^-\-<X) 


5929 
4665600■ 


We continue the above correction process to successively determine the correction function 
MCfc( 3 :) until some k* you want. On one hand, to find the related coefficients, we often use 
an appropriate symbolic computations software because it’s huge of computations. On the 
other hand, the exact expressions at each occurrence also need lot of space. Hence in this paper 
we omit many related details. For interesting readers, see our previous paper [SIE]. In fact, we 
can prove that for 0 < fc < 3 


(3.8) 

where 


Kq = 


Ki = 


K2 = 


l<'3 = 


MCfc(x) = K 


i—0 X + A,- 


1 

72’ 

5929 

32400’ 

76899172249 

248039857296’ 
786873417270631211749921 


Ao — 


31 


Ai — 


90’ 
481937 


A 2 — 


3735270’ 
7745462509019287 


A 3 — 


19149278075101482’ 
2098335745817751685364201067279071 


851541507731717527392144’ " 30311088872486921466334781589254970 

By Lemma 1 again, we get for some constant Ck ^ 0 

(3.9) lim x‘^^^^Ek{x) = Ck, (/c = 0,1,2, 3), 

x—>-+oo 


i.e. n{Ei^{x)) = 2A; + 4 for A: = 0,1, 2 , 3. Thus we obtain more accurate approximation formulas: 


(3.10) T{x + l) = V^{^ + ^ + MCfc(a:) • (^1 + 0{x (^^+4))^ ^ 3 ;; ^ _|_oo. 

It should be noted that if we rewrite MCfc(x) in the form of where P, Q are polynomials 

with r = k and s = A: + 1, theoretically at least, for a large x the above formula may reduce or 
eliminate numerically computations compared with the previous results, see e.g. [9l[T2]. This is 
the main advantage of the multiple-correction method. 

The following theorem tells us how to obtain sharp inequalities. 

Theorem 1. Let MCfc(x) be defined as (13.81) . Let x>l, then we have for k = 0,2, 


(3.11) r(x + 1) > AlCfc(a;), 

and for k = 1,3, 

(3.12) T{x-\-l) < + ^ + MCfc(x). 
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Proof. We let fk{x) = Ek{x) — Ek{x + 1). Clearly if liuix^^oo Ek{x) = 0, then Ek{x) = 
+ j). This transformation plays an important role in this paper (essentially, it is 
a difference method). Hence, in order to prove inequality Ek{x) > 0 (or Ek{x) < 0), it suffices 
to show that the equality fk{x) > 0 (or fk{x) < 0 ) holds under the condition lima,^+oo Ek{x) = 0 . 
By the Stirling’s formula (1.2), we can show that the condition limj,_>._|_oo Ek{x) = 0 always holds. 
In what follows, we will apply this condition many times. 

By using Mathematica software, we may prove that for x > 1 


foix) 

fiix) 

fsi^) 


_^i( 8 ;x)_ 

x(l + x)2(31 + 90x)2(121 + 90x) 2(77 + 552x + 1080x2)2(1709 + 2712x + 1080x2)2 
^2(13;x)(x - 1) + 1463••• 9447 
~x(l + x) 2)(1359251 + 2829648x + 5976432x2)2^3(16; x) ’ 

^ 4 ( 20 ; x) 

x(l + x) 2 T 5 ( 28 ;x) 

^6(25; x)(x - 1) + 17135 • • • 66999 ^ 

x(l + x) 2T7(36; x) ^ 


> 0 , 


We only give the proof of inequalities in case k = 3, other may be proved similarly. In this case, 
we see that for x > 1 the inequality (I3.12p is equivalent to £' 3 (x) < 0. As lim 3 ;^+oo £■ 3 ( 2 ;) = 0, it 
suffices to prove that fsix) < 0 for x > 1. Since /^(x) is strictly decreasing, but lim 2 ;_>.+oo fsix) = 
0, so f^{x) > 0. Thus fsix) is strictly increasing with liuix^+oc fsix) = 0, so fsix) < 0. This 
completes the proof of Theorem 1. □ 

By the multiple-correction method, we also find another kind of inequalities. 

Theorem 2. Let the k-th correction function MCfc(x) be defined by 


MCo(x) 


MCfc(x) 


KO 


{x + 1^)2 + Ao 


Kq 


k 

K 


(x + T ^oT ^ 4“ 


{k > 1 ), 


where 


Ko = 


Kl = 


K2 = 


1 

“144’ 

4394 

637875’ 
7894414898425 


Ai — 


119793516544 ’ 

1897560849252106177858465792 
77174813342532578267347147395 ^ 

Then we have 


= 


A 3 — 


4007 

° “ 21600’ 
130311599 

15575040 ’ 
265702682899837009577 

34427631789478287360 ’ 
30320380455616293004898928163131563244811979 

6134364315672065325746652708240298034227200 ’ 


A2 — — 


(3.13) 


r(x +1) < (— 

V e 


X + - (1 + MCo(x)), X > 13, 
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(3.14) 


T{x + 1) < yx + ^ (1 + MC2 (x)) , X > 6, 

and for k = 1,3, 

(3.15) r(x + 1) > \/^ ^—^ y^x + - (1 + MCfc(x)), X > 1. 

Proof. Since the proof of Theorem 2 is very similar to that of Theorem 1, here we only give the 
outline of the proof. First, let the relative error function Ek{x) be dehned by 

(3.16) T{x + l) = V^{^ ^x + i (1 + MCfc(x)) exp(F;fc(x)). 

Hence 

(3.17) Et(x) - + 1) = -1 + iln (l + i) + In 

By making use of Mathematica software and Lemma 1, we can prove 

(3.18) /r(£’fc(x)) = 2A: + 5, (A: = 0,1, 2, 3). 

Next we let gk{x) = Ek{x) — Ek{x + 1). By using Mathematica software, it isn’t difficult to check 
that 


9q{x) 

g'iix) 

92{x) 

g'iix) 


4'i(14; x)(x - 13) + 29707 • • • 81369 ^ 
x(l + x)2(l + 6x)2(7 + 6 x) 24'2(16; x) 
^3(20; x)(x - 1) + 13798 • • • 89479 
x(l + x)2(l + 6x)2(7 + 6 x)2'I'4(24; x) 
^5(26; x)(x - 6) + 97250 • • • 34321 ^ 

x(l + x)2(l + 6x)2(7 + 6 x) 24'6(32; x) 
^7(32; x)(x - 1) + 836559 • • • 37479 
x(l + x)2(l + 6x)2(7 + 6 x)2'I'8(40; x) 


Lastly, just as the proof of Theoreml, Theorem 2 follows from the above inequalities readily. □ 


4 Ramanujan-type inequalities 

Theorem 3. Let the k-th correction function MCfc(x) be defined as 


(4.1) 


MCfc(x) = K 


j=o X + 6,- ’ 
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where 


11 

““ - 

459733 

dl — -, 

711480 

49600874140433 

“ 101450127018720 ’ 

_ 169085305336152527131511003963 ^ 

“ 101221579151797375403194730976 ’ ^ 

Let X > 1, then for k = 0,2, 


1455925 

^ ~ 70798882 ’ 

_ 10259108965771635091 

^ “ 19545564575317443762 ’ 
6141448535908002711219920016488834171 

203275987838924050801436670299517447102' 


(4.2) 


r(x + l)<v^^^^ ^8x^ + 4x^ + X + ^ + MCfc(x)^ , 


and for k = 1,3, 


(4.3) 


r(x +1) > (- 

\ e 


X 


8x^ + 4x^ + X + ■— + MCfc(x^ 

oU 


i 

6 


Proof. We define the relative error function Ek{x) by the relation 

(4.4) T{x+ 1) = y/f: f8x^ + 4x^ + X + ^ + MCfc(x)^ exp(£’fc(x)). 


Thus 


(4.5) 


Ek{x) - Ek{x + l) = -l + xln( 1 + - 

X 


1 8(x + 1) + 4(x + 1) + (x + 1) + -^ + MCfc(x + 1) 

H— in- 

6 


8 x^ + 4x^ + X + ^ + MCfc(x) 
By using Mathematica software and Lemma 1, we can check 


(4.6) 


^i{Ek{x)) = 2k + 0, {k = 0,1,2,3). 


We let Uk{x) = Ek{x) — Ek[x + 1). By making use of Mathematica software again, we can prove 


U'^{x) 

U'{{x) 

U'f{x) 

U'i{x) 


_Ti(13;x)(x - 1) + 416838558509297754261614731715717_ 

3x(l + x)2(79 + 154x)2(233 + 154x)2 (T2i(3; x)(x - 1) + 363565)^ ^ 22 ( 4 ; x) 
T3(19 ;x)(x - 1) + 85653-••25001 

x(l + x) 2 'I' 4 ( 28 ; x) ^ ’ 

4^5(25; x)(x - 1) + 32968 • • • 13479 

x(l + x) 24'6(36; x) ^ ’ 

4'7(31; x)(x - 1) + 17145 • • • 57723 

3x(l + x) 24'8(44; x) ^ 


Similar to the proof of Theorem 1, we can get the desired assertions from the above inequalities. 

□ 













Theorem 4. Let the first-correction function MC^(x) be defined by 


(4.7) 

where 


MCt(x) = 


Ko 


+ 


Kl 


x + Ao + Aiox^ + Aiix + A 


12 


11 

459733 

Kl = 


Ao — 


79 


All — 


15523200 ’ 
717183502490887 


Am — 


154’ 
71181889 


520777318696096 
Then for x >1, the following inequality holds true 


Ai2 — 


70798882’ 
1118629052995381153799 

1958878792277282473920’ 


(4.8) r(x + l)<v^^^^ ^8x^ + 4x^ + X + ^ + MCi(x) 

Proof. Let the first-correction error function E*(x) be defined by 


(4.9) 

Hence 


r(x-7 1) = 0r f 8x^-7 4x^-7 X-h ^-h MC*(x) j exp(L;j'(x)) 


(4.10) .Ej'(x)-L;i*(x + l)=-l + xln 1 + - 


1 ^ 8(x -|- 1)^ -|- 4(x -|- 1)^ -|- (x -|- 1) ^ MC^(x -|- 1) 

6 Sx^ -|- 4x2 -|- X ^ MCi(x) 

By using Mathematica software and Lemma 1, we have 
(4.11) h{El{x)) = 10. 

Now we let V{x) = L^i(x) — Lij['(x -|- 1). By using Mathematica again, we have 

4-1 (33; x)(x - 1) -h 96057 • • • 27429 


(4.12) H"(x) = - 


3x (4-2(3; x))^ 4-3(12; x) ( 4 - 4 ( 6 ; x)(x - 1) -7 2169 • • • 3461)^ 4-5(14; x) 


< 0. 


By the same approach as the proof of Theorem 1, the inequality (j4.8p follows from the (j4.12h . □ 

Remark 2. It is an interesting question whether our method may be used to obtain some sharp 
bounds for the ratio of the gamma functions, see e.g. [IIl[l9l|20l[2l]. 
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